Renormalization of the tunnel splitting in a rotating nanomagnet 
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We study spin tunneling in a magnetic nanoparticle with biaxial anisotropy that is free to rotate 
about its anisotropy axis. Exact instanton of the coupled equations of motion is found that connects 
degenerate classical energy minima. We show that mechanical freedom of the particle renormalizes 
magnetic anisotropy and increases the tunnel splitting. 
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Macroscopic dynamics of a fixed-length magnetic mo- 
ment, M, of a single-domain ferromagnetic particle is 
described by the Landau-Lifshitz equation 1,2 . When dis- 
sipation (which is usually weak) is neglected this equation 
reads 
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where £ is the classical magnetic energy of the particle 
that depends on the orientation of M. It was shown long 
ago^ that Eq. ([T]), besides the real-time solutions, also 
possesses imaginary-time solutions - instantons - that 
describe macroscopic quantum tunneling of M between 
classically degenerate energy minima (see also books: 
Refs. and l5|) . In early experiments on spin tunneling^ 
single-domain magnetic particles were always frozen in 
a solid matrix so that their physical position and orien- 
tation were fixed and only rotation of the magnetic mo- 
ment was allowed. Later, beams of small magnetic clus- 
ters were investigated^— and more recently free mag- 
netic nanoparticles confined within solid nanocavities 
have been studied^. Experimentalists have also worked 
with molecular nanomagnets deposited on surfaces^— 
or carbon nanotubes 17 , as well as with single magnetic 
molecules bridged between metallic electrodes^—. In 
such experiments the particles retain some mechanical 
freedom. This inspired recent theoretical work on quan- 
tum mechanics of rotating magnets^—. In this paper 
we obtain exact magnetic instanton for a single-domain 
particle that is free to rotate about its anisotropy axis. 

General analytical solution for the rotational quantum 
levels of a rigid body does not existed. Spin degree of 
freedom makes this problem even less tractable. How- 
ever, as was recently demonstrated in Ref. the exact 
eigenstates and exact energy levels can be obtained ana- 
lytically for a nanomagnet that, due to a large magnetic 
anisotropy, can be described as a two-state spin system 
and is free to rotate about its magnetic anisotropy axis. 
Such a system at rest is described by the Hamiltonian 
T~i a = — (A/2)a x , where a x is a Pauli matrix and A is 
the tunnel splitting of spin-up and spin-down states, with 
z being the anisotropy axis. It arises from the terms in 
the full Hamiltonian that do not commute with the z- 
component of the total spin S. Rotation of the particle 
by the angle tp about the Z-&xis transforms the two-state 
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FIG. 1. Color online: Spin and angular momentum tunneling 
in a rotating magnetic nanoparticle. 



Hamiltonian into 2 ^ 
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(S is dimensionless). Exact eigenvalues of the full Hamil- 
tonian of a nanomagnet rotating about its anisotropy 
axis, H = (HL Z ) 2 /(2I) + H' a (where L z is the z- 
component of the dimensionless mechanical angular mo- 
mentum and I is the moment of inertia), were obtained 
in Ref. l23l . where it was shown that parameter 



a = 2(HS) 2 /(IA) 
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determines low energy states of the particle. At a < 
oti = [l — l/(25) 2 ] the ground state and the first ex- 
cited state are respectively symmetric and antisymmetric 
superpositions ofJ = L + S = states shown in Fig. [TJ 
with energies E± — ^f - ± —■ 

Derivation of these results was based upon the assump- 
tion that the parameter A is the same for a stationary 
magnetic particle and for a particle that is free to rotate. 
Instanton method allows one to test this assumption. Be- 
low we find the exact instanton solution of the equations 
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fined with respect to the particle-frame axes, 
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Here /xo is the magnetic permeability of vacuum, V is the 
volume of the particle and Mq = M/V represents the 
magnetization which is a constant of the ferromagnetic 
material. Anisotropy constants K± > and K\\ > have 
been redefined to show the explicit proportionality of the 
anisotropy energy to the volume. They are dimensionless 
numbers, typically of order unity. 

The rotational kinetic energy of the particle is £r = 
\l<p 2 ■ The Lagrangian of our system consists of the triv- 
ial kinetic Lagrangian for the variable ip, 



FIG. 2. Geometry of particle and laboratory frames. 



Cl = KL ■ ip 

and the magnetic Lagrangian 



(5) 



of motion describing the dynamics of the magnetic mo- 
ment and the rotation of the particle. It shows that me- 
chanical freedom does renormalize the tunnel splitting A. 
However, this renormalization is small unless A is very 
large and a is close to a\. 

Consider a high-spin magnetic particle with biaxial 
anisotropy that is free to rotate about its easy axis. The 
initial state of the particle is such that its total angular 
momentum is zero, J = S + L = 0. In other words, the 
total spin (magnetization) vector points along the easy 
axis and the particle rotates about this axis such that 
these angular momenta are equal in magnitude and oppo- 
site in direction, see Fig. [TJ The exchange interaction be- 
tween individual spins is strong, so the magnitude of the 
total spin of the particle is a constant. The magnetic en- 
ergy will be expressed below in terms of M which is pro- 
portional to the spin, hS = M/7. Here 7 = — e/2mc < 
is the electron gyromagnetic ratio. The orbital angular 
momentum is associated with the rotational motion of 
the particle itself, KL = lip, where I is the particle's 
moment of inertia and <p is its angular velocity. 

We define coordinate systems of the lab frame (x, y, z) 
and particle frame (X, Y, Z) as shown in Fig. [2] In the 
particle frame the x-axis is along the easy axis in the xy 
easy plane, and the z-axis is the hard axis. The lab frame 
is centered at the same origin such that the X-axis of the 
lab frame coincides with the x-axis of the particle frame. 
The particle is free to rotate about this axis. At some 
initial time t = we choose the two coordinate frames to 
coincide. The angle of rotation of the x- and z-axes with 
respect to the X- and Z-ax.es is <p(t). Notice the change 
of the easy axis as compared to the choice of Ref. |23| . 
which is dictated by the mathematics of the problem. 

The anisotropy energy is naturally defined in the par- 
ticle frame: £a = k±M 2 — k\\M%. It can be written in 
terms of spherical polar coordinates {9, 0) which are de- 
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for the variables 9 and <fi. The first term in Eq. ([6]) fol- 
lows from the fact that HS Z = hS cos 9 is the generalized 
momentum for the coordinate </>. The second term is the 
effective magnetic energy in the rotating frame, 



£' A =£ A (9,<j))-hS-cp. 
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The last term in this equation is related to the fact that 
in the particle frame the rotation is equivalent to the 
magnetic field B = 93/7. 

The total Lagrangian of the particle is a sum of Cl 
and Cs- 



C = H(L + S)-(p- 
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The first term reflects the fact that in the presence of a 
spin the generator of rotations is J = L + S. The explicit 
form of the total Lagrangian in terms of the generalized 
coordinates 9, 4>, and (p is 
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The equations of motion are Euler-Lagrange equations 
for 9, <f>, and <p: 
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where we have introduced dimensionless time t — 
jfioMot and (p = dtp/dt. 

Note that the equations of motion for <f> and 9 can 
also be obtained from the Landau-Lifshitz equation with 
£ = £' A . Indeed, the equations for <fi and 9 that follow 
from Eq. (see, e.g., Ref.0), 
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are identical to the equations (fTU| and (fTTj). The third 
equation of motion, Eq. \12\ is the conservation of the 
total angular momentum: 
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At Jx = it is equivalent to the constraint: 
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With account of this constraint the equations of motion 
for d> and 9 become 



= (K ± + K cos 2 0)cos0 
dt " 

d(cos9) . 2 
= — = — Ku cos © sin « sin I 
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We see that for J = the effect of rotations reduces 
to the renormalization of the easy axis anisotropy ATy. 
This is easy to understand from the following considera- 
tion. In a state with J x = 0, equilibrium vectors S and 
L look in the opposite directions along the s-axis. If 
S deviates from the x-axis, S x decreases and so should 
L x to preserve the condition J x — 0. The decrease of 
L x corresponds to the decrease of the rotational energy, 
(HL X ) 2 /(2I), mandated by J x — 0. Thus, effectively, the 
magnetic anisotropy energy associated with the deviation 
of S from the easy axis becomes smaller when mechanical 
rotation is allowed. 

We should now look for solutions of equations (fl!))) and 
(fTf). We hrst notice that 

£ = ln M^V[(K ± + K'u cos 2 4>) cos 2 9 + K[ sin 2 0] (19) 



is the integral of motion. This is easy to see by differen- 
tiating this equation on time and substituting in the re- 
sulting equation the time derivatives <f> and 9 from equa- 
tions (|16|) and (IT7|) . Not surprisingly, up to a constant, 
Eq. ([T9"]) equals the total energy of the particle in the 
laboratory frame, £ — £a + ^Itp 2 , with account of the 
constraint (|T5|) . Eq. (fTTj)) gives 
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FIG. 3. Dependence of on r for clockwise and counterclock- 
wise underbarrier rotations of M. 



This allows one to express £ in terms of the angle 
its time derivative: 
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Since this expression is positively defined, the classical 
energy minima occur at £ = 0. They correspond to 
the stationary magnetization pointing in either direction 
along the easy axis, i.e., <f> = 0> ?r with cos 9 = in accor- 
dance with Eq. ((211)) . 

Equation £ = has no real-time solutions for (f> 
that connect the two degenerate classical energy minima. 
However, in imaginary time, f = it, equation £ = is 
equivalent to 
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Such equation has instanton solutions that connect the 
classical energy minima: 



(t) = iarccos < — 
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The T-dependence of 9 is given by Eq. (|2U|) . and the r- 
dependence of tp is given by Eq. (1151) . Here 



\ = K{/K U tj = K'^K'u + K ± ) . (24) 

The positive and negative signs correspond to the two 
possible trajectories, which are counterclockwise and 
clockwise rotations of the magnetization from <p = at 
t = — oo to 4> — i 71 " at r = +oo, respectively, see Fig- EI 

The tunnel splitting has the form A = Ae B , where A is 
of the order of quantized oscillations near the minimum 
of the potential well and 
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is the WKB exponent. Substituting here C of Eq. ([8} at 
J = 0, one obtains for the instanton trajectory 
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To see the effect of the mechanical freedom of the parti- 
cle on spin tunneling we define a dimensionless parameter 
a' = Kr/Ku. In a microscopic theory the easy-axis crys- 
tal field is presented as —DS 2 . The connection between 
D and the parameter K\\ of the macroscopic theory is 2 



(27) 



where s is spin per unit cell of the crystal and Vq is the 
volume of the unit cell. (Singularity at s = 1/2 reflects 
the fact that single-ion magnetic anisotropy does not ex- 
ist for spin 1/2). The total spin of a ferromagnetic par- 
ticle can be presented as S = s(V/Vq). Consequently, 
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where we have used Eqs. ([3]) and (|18p . Renormalization 
of the easy-axis anisotropy by rotations can be presented 
in the form K!, = K»(l — e), where 
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and Ei — (2S— l)D is the energy of the first excited spin 
state at A = 0. The low energy limit that we are studying 
corresponds to A < i?i and a < ai = [l — 1/(2S) 2 ] , 
see Ref. [23|. In this limit e is small. Consider, e.g., the 
case of large S and large A (small tunneling rate). Ac- 
cording to Eq. ([2"tJ| in this case B = — 51n(4A) so that 
A oc exp [— 51n(4A)]. It is easy to see from this expres- 
sion that mechanical rotation renormalizes A by a factor 
exp (eS) . Normally it would not be large compared to 
one. However, since small e in the exponent is multiplied 
by a large S, it is not out of question that at sufficiently 
large A a slight increase of the tunnel splitting would be 
observable in spin clusters that are free to rotate. 
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